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Galois conjugation relates unitary conformal field theories (CFTs) and topological quantum field theories 
(TQFTs) to their non-unitary counterparts. Here we investigate Galois conjugates of quantum double models, 
such as the Levin- Wen model. While these Galois conjugated Hamiltonians are typically non-Hermitian, we find 
that their ground state wave functions still obey a generalized version of the usual code property (local operators 
do not act on the ground state manifold) and hence enjoy a generalized topological protection. The key question 
addressed in this paper is whether such non-unitary topological phases can also appear as the ground states of 
Hermitian Hamiltonians. Specific attempts at constructing Hermitian Hamiltonians with these ground states 
lead to a loss of the code property and topological protection of the degenerate ground states. Beyond this we 
rigorously prove that no local change of basis l |IV.5[ > can transform the ground states of the Galois conjugated 
doubled Fibonacci theory into the ground states of a topological model whose Hermitian Hamiltonian satisfies 
Lieb-Robinson bounds. These include all gapped local or quasi-local Hamiltonians. A similar statement holds 
for many other non-unitary TQFTs. One consequence is that the "Gaffnian" wave function cannot be the ground 
state of a gapped fractional quantum Hall state. 

PACS numbers: 05.30.Pr, 73.43.-f 



I. INTRODUCTION 

Galois conjugation, by definition, replaces a root of a poly- 
nomial by another one with identical algebraic properties. For 
example, i and — i are Galois conjugate (consider z 2 + 1 = 0) 

as are = and - \ = (consider z 2 - z - 1 = 0), 

as well as \/2, v^e 2 ™/ 3 , and S^e" 2 ™/ 3 (consider z 3 - 2 = 
0). In physics Galois conjugation can be used to convert non- 
unitary conformal field theories (CFTs) to unitary ones, and 
vice versa. One famous example is the non-unitary Yang-Lee 
CFT, which is Galois conjugate to the Fibonacci CFT (G?2)i, 
the even (or integer-spin) subset of su(2) 3 . 

In statistical mechanics non-unitary conformal field theo- 
ries have a venerable history.^ However, it has remained less 
clear if there exist physical situations in which non-unitary 
models can provide a useful description of the low energy 
physics of a quantum mechanical system - after all, Galois 
conjugation typically destroys the Hermitian property of the 
Hamiltonian. Some non-Hermitian Hamiltonians, which sur- 
prisingly have totally real spectrum, have been found to arise 
in the study of PT-invariant one-particle systems 3 and in 
some Galois conjugate many-body systems 4 and might be 
seen to open the door a crack to the physical use of such 
models. Another situation, which has recently attracted some 
interest, is the question whether non-unitary models can de- 
scribe ID edge states of certain 2D bulk states (the edge holo- 
graphic for the bulk). In particular, there is currently a discus- 
sion on whether or not the "Gaffnian" wave function could be 
the ground state for a gapped fractional quantum Hall (FQH) 
state albeit with a non-unitary "Yang-Lee" CFT describing its 
edge P™^ We conclude that this is not possible, further restrict- 
ing the possible scope of non-unitary models in quantum me- 
chanics. 

We reach this conclusion quite indirectly. Our main thrust 
is the investigation of Galois conjugation in the simplest non- 



Abelian Levin- Wen model. This model, which is also called 
"DFib", is a topological quantum field theory (TQFT) whose 
states are string-nets on a surface labeled by either a triv- 
ial or "Fibonacci" anyon. From this starting point, we give 
a rigorous argument that the "Gaffnian" ground state cannot 
be locally conjugated to the ground state of any topological 
phase, within a Hermitian model satisfying Lieb-Robinson 
(LR) bounds 9 (which includes but is not limited to gapped 
local and quasi-local Hamiltonians). 

Lieb-Robinson bounds are a technical tool for local lattice 
models. In relativistically invariant field theories, the speed of 
light is a strict upper bound to the velocity of propagation. In 
lattice theories, the LR bounds provide a similar upper bound 
by a velocity called the LR velocity, but in contrast to the rel- 
ativistic case there can be some exponentially small "leakage" 
outside the light-cone in the lattice case. The Lieb-Robinson 
bounds are a way of bounding the leakage outside the light- 
cone. The LR velocity is set by microscopic details of the 
Hamiltonian, such as the interaction strength and range. Com- 
bining the LR bounds with the spectral gap enables us to prove 
locality of various correlation and response functions. We will 
call a Hamiltonian a Lieb-Robinson Hamiltonian if it satisfies 
LR bounds. 

We work primarily with a single example, but it should be 
clear that the concept of Galois conjugation can be widely ap- 
plied to TQFTs. The essential idea is to retain the particle 
types and fusion rules of a unitary theory but when one comes 
to writing down the algebraic form of the F-matrices (also 
called 6j symbols), the entries are now Galois conjugated. A 
slight complication, which is actually an asset, is that writing 
an F-matrix requires a gauge choice and the most convenient 
choice may differ before and after Galois conjugation. 

Our method is not restricted to Galois conjugated DFib e 

and its factors Fib e and Fib e , but can be generalized to in- 
finitely many non-unitary TQFTs, showing that they will not 
arise as low energy models for a gapped 2D quantum mechan- 



ical system with topological order. 

The 2D quantum mechanical systems which can be de- 
scribed by any type of TQFTs are known as topological 
phases. Although this concept is widely noted in the con- 
densed matter physics literature, our introduction is not com- 
plete without providing a definition. Many authors focus on 
properties (e.g. existence of anyonic excitations), but we pre- 
fer to give a more fundamental definition since it is this defi- 
nition which figures into our proof in Sec. |IV] We say that a 
LR Hamiltonian describes a topological phase (or a phase is 
topologically ordered) if and only if its ground state manifold 
G satisfies the following "code" property with respect to all 
spacially local operators L: the composition 



G^H 



H 



-> G 



(1) 



is a multiplication by some scalar s(L) (possibly s(L) = 0). 
L local means L acts only on sites in a sufficiently small ra- 
dius. This definition first appeared in Ref. [Tol (see definition 
3.6 there), conceptualizing the earlier formulation of Ref. [TT1 
Recently Bravyi et a/M^lhave called this axiom "TQO-1" and 
advocated an additional requirement that they called "TQO- 
2", which enforces a consistency between local and global 
ground states. While technically necessary, we know of no 
realistic case where the second axiom would be required and 
so have not included it in the framework of this paperfSl 



H. LEVIN-WEN MODEL AND ITS GALOIS CONJUGATES 
A. The Levin- Wen model 

Topological quantum field theories are highly constrained 
mathematical construct^GH designed to capture the low en- 
ergy physics of topologically ordered systems. Chern-Simons 
theory-^ generates most of the known examples; the simplest 
of these, all chiral, being based on a Lie group and level k, 
Gfc. Starting from a set of particles and fusion rules, there is a 
standard construction - called the "quantum double" or "Drin- 
feld center" - which produces an achiral TQFT Such quan- 
tum doubles were introduced in the physics literature by Levin 
and Wen 8 in the form of "string-net" Hamiltonians. If, for in- 
stance, we take the particles and fusion rules from the chiral 
Fib TQFT, see Eq.(|2]i below, and use these to label string-nets 
on surfaces, a "larger" TQFT DFib = Fib <E> Fib (with more 
particle types) is obtained. 

The Levin-Wen model thus is a microscopic spin Hamilto- 
nian implementing doubled topological theories. Originally, 
it was defined 8 -! on a honeycomb lattice, but its extension to 
any trivalent graph is straight-forward. Given a lattice graph 
and an anyonic theory, the model's Hilbert space is spanned 
by all labelings of graph edges with the theory's particle types 
which are consistent with a set of constraints given by the the- 
ory, the so-called fusion rules. As a simple example we first 
consider the Fibonacci theory Fib, where there are only two 
particle types, namely a trivial particle 1 and the Fibonacci 
anyon r. Two particles can combine according to the fusion 




FIG. 1. Edge labeling for a plaquette of the honeycomb lattice. 



rules 



1x1 = 1 



1 X T 



TXT 



1 + T 



(2) 



In the Levin- Wen model implementing the doubled Fibonacci 
theory DFib, this amounts to the constraint that of the three 
edges meeting in any single vertex never only one can carry a 
t label. This Hilbert space can either be understood as that of 
an anyonic quantum liquid enclosing the lattice links^ or al- 
ternatively as the the ground states of a spin model (by identi- 
fying particle types with spin directions) with a peculiar three- 
spin interaction enforcing the vertex constraint. 



Within these states, the Hamiltonian 



H LW = J, 



plaquettes p 



(3) 



is a projector onto the r-flux state of a plaquette p thus fa- 
voring the trivial flux cf>(p) = 1 through each plaquette. The 
action of this operator on an element of the basis where the 
edges belonging to plaquette p carry labels a, . . . , a, . . . , f 
as displayed in Fig.[T]results in a superposition of states where 
the inner edges of the plaquette carry new labels a', . . . , £' 
whereas all other edges remain unchanged. Any of the labels 
takes one of the values {1, t}. The matrix elements between 
these basis states read explicitly (see Refs.|8land[T8lfor a de- 
tailed derivation) 



k-); (*-): 

s 



where d s denotes the quantum dimension of particle type s, 
i.e. e?i = 1 and d T = 4> = (1 + \/5)/2, the golden ratio and 
D the total quantum dimension, D = \J d\ + d\ = \J2 + (j> 
for Fibonacci anyons. For different plaquette geometries this 
operator has an analogous form with one F-symbol for each 
edge of the plaquette. 

The F-symbol, which can be thought of as a generalized 6j 
symbol, describes local basis transformations in a fusion tree 
as shown in Fig. [2] and is a defining property of the anyonic 
theory. For Fibonacci anyons, this transformation is trivial 
except for the case when all four outer legs of the subgraph 
that is to be transformed carry the r label. Then we have the 
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unitary 2x2 matrix 



Fl 



0- 1 c^ 1 ' 2 



(5) 



where any choice of the gauge c satisfies the pentagon equa- 
tions for the Fibonacci fusion rules. Choosing the gauge =1 
gives the usual unitary F-matrix for the Fibonacci theory, 
which we refer to as the symmetric normalization. From an al- 
gebraic point of view the natural gauge choice is =0 5//2 which 
leads to 



Fl 



0-1 4> + i 

20-3 1-0 



(6) 



where no square roots of appear. We refer to this choice 
as the algebraic normalization and define A = c/0 5 / 2 . We 
remark that both normalizations give the same spectra for our 
models, since the corresponding Hamiltonians are conjugate 
by a diagonal fugacity change matrix. 

The Levin-Wen model can can be solved exactly since 
all the plaquette terms commute. 8 ' As a sum of projectors it 
counts the number of plaquettes penetrated by a nontrivial 
T-flux and the spectrum hence consists of states at all non- 
negative integer multiples of J p , corresponding to the number 
of nontrivial plaquette fluxes. 

The ground states of the model correspond to all states 
with no plaquette fluxes, corresponding to the ground states 
of the topological liquid on a doubled surface around the lat- 
tice. With periodic boundary conditions in both directions this 
surface is a doubled torus with four degenerate ground states. 



B. The doubled Yang-Lee model 

Now we turn to a theory of non-unitary non-abelian anyons 
which are closely related to the Fibonacci ones by Galois con- 
jugation. We start by noting that Fib is only one particular 
sub-theory out of a discrete set of su(2)k (for finite k) any- 
onic theories, specifically the integral spin half of the unitary 
su(2)3 theory. The su(2)k theories are certain deformations 
of SU(2) characterized by the truncation level k, which de- 
fines the particle types in the theory, and additionally a de- 
formation parameter q determining the precise values of the 
F-symbol. For fc = 3 two sub-theories characterized by the 
roots of unity t = e 27 ™/ 5 and t' = e 47 ™/ 5 are Galois conju- 
gates of each other (see Fig. [3}. 

Now, the former value for q produces the Fibonacci theory 
as described above whereas the latter leads to the non-unitary 
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FIG. 2. The F -symbol. 
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FIG. 3. (color online) The i-deformation parameters of Fibonacci 
and Yang-Lee anyons correspond to different primitive roots of unity. 



F matrix 



FZ 



d-1 d+l 
2d - 3 1 - d 



in the algebraic normalization and 



-0 icj) 1 ' 2 

i(t) 1/2 



(7) 



(8) 



in the symmetric normalization. These are just the F-matrices 
of the DFib theory with = — t 1 ' 2 — t~ x l 2 replaced by d = 
— 1/0 = — t' 1 / 2 — t'~ x l 2 . Here we choose the fourth roots 
f i/4 = j e Wio and t n/A = ie -37ri/i0 ) which w01 be needed 

to specify the Galois conjugation of the full theory below. 

We remark on the choice of algebraic normalization for the 
F matrix. For the Yang-Lee theory no choice of c would make 
F unitary — a manifestation of the non-unitarity. While there 
would be no topological invariant positive definite Hermi- 
tian products on all ground state manifolds in Yang-Lee the- 
ory, there is always a topological invariant Hermitian product 
with possibly mixed signatures. The topological invariant in- 
ner product is Hermitian because the partition function under 
time reversal in a (2+l)-topological theory is Hermitian conju- 
gated. The (1, 2)-entry of the above F-matrix is the theta sym- 
bol (the norm of a fusion basis in a fusion space) multiplied 
by dr 2 , hence it should be a real number. The above choice 
of F-matrix for Fibonacci case is pleasant when we work on 
number theory related problems. In particular, one notices 
that in this case we obtain the Galois conjugate by replacing 
all occurrences of the golden ratio by d = —1/0, which is 
the second solution of the quadratic equation x 2 = 1 + x. 

As Galois conjugation does not change the theory's al- 
gebraic structure the doubled Yang-Lee (DYL) Levin-Wen 
model using the F matrix of Eq. ([SJ, can be solved in exactly 
the same way as its DFib counterpart. In particular, it has 
exactly the same spectrum whose eigenvalues count the num- 
ber of plaquettes penetrated by a non-trivial flux and the same 
ground state degeneracies. The DYL model also retains the 
topological protection of the ground state degeneracy against 
local perturbations. 
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III. HERMITIAN MODEL FROM NON-UNITARY 
THEORY 

A. Constructing Hermitian models 

While the non-Hermitian DYL model features a general- 
ized stable topological phase and a generalized code property, 
discussed in more detail below, an immediately arising ques- 
tion is whether this phase can also be realized in a Hermitian 
model. There are multiple ways to obtain a Hermitian model 
that has the same ground states as the non-Hermitian parent 
model. However, as we will see in the following the question 
whether the topological nature of the ground state remains is 
a more subtle one. 

The simplest Hermitian model H is obtained by squar- 
ing the non-Hermitian parent Hamiltonian H. This model 
has the same right ground-state eigenvectors as the original 
model. Alternatively, HH^ has the same left ground-state 
eigenvectors. The simplicity of this approach comes at the 
cost of a Hamiltonian which is highly non-local. To avoid 
non-local terms, we can take an alternative route and individ- 
ually square each plaquette term of H p — <5°^ T , arriving 

at the Hamiltonian J2 P Hj,H p or J2 P H p H^. Since each pla- 
quette term annihilates the ground state, squaring them in this 
way also annihilate the (right/left) ground state eigenvectors. 
Finally, we can replace the non-Hermitian plaquette operator 
H p with a projector onto the complement of the operator's 
kernel. More specifically, we diagonalize the plaquette oper- 



ator and use its orthogonalized right eigenvectors 
longing to the eigenvalue to define a projector 



(r) 



be- 



Independent of the way we derive a Hermitian model from 
the parent DYL model, we find that the code property is lost 
for the Hermitian models: when constructing a Hermitian 
model, one inevitably has to decide wether to preserve left 
or right ground states. The code property for the Hermitian 
model would require expectation values of local operators of 
the form 



(r) 



(r) 



and 



(0 



(12) 



to again be multiples of the identity. In general, this usual 
code property will not be satisfied, as one can see, for exam- 
ple, by calculating the matrix elements of a local observable 
such as a string tension. Perturbing any Hermitian Hamilto- 
nian which has the (right or left) DYL ground states with an 
arbitrary small string tension will hence immediately lead to 
a splitting of the ground-state degeneracy, as we will discuss 
below. 



C. Absence of topological order 

In this section we probe whether topological order survives 
the construction of a Hermitian model by numerically diago- 
nalizing the models on different lattice geometries, the hon- 
eycomb lattice of the original Levein-Wen constructiorP and 
the two-leg ladder geometry of Ref. 18 We diagonalized sys- 
tems with up to 24 edges using a dense eigenvalue solver and 
employed iterative schemes for systems with up to 39 edges: 
the Lanczos algorithm for Hermitian models and an implicitly 
restarted Arnoldi method for non-Hermitian models. 



1 



E 



(r) 



(r) 



(9) 



1. Honeycomb model 



The sum of these projectors is then used to define the Hermi- 
tian Hamiltonian 



H 



herm 



(10) 



It turns out that all three approaches result in the same qualita- 
tive behavior - a loss of the code property and the associated 
stable topological order - and we will limit our discussion to 
the last approach. 



Our results on the honeycomb lattice show a clear distinc- 
tion between the DFib and DYL models on the one hand and 
the Hermitian model iJ horm derived from the DYL model 
on the other hand. While all models feature four degenerate 
ground states, the former two are gapped, whereas the latter 
one turns out to be gapless in the thermodynamic limit; see the 
finite-size extrapolation in Fig.|4^). Furthermore, the ground- 
state degeneracy is easily lifted by a local perturbation, such 
as a string tension - in contrast to the stability of the topolog- 
ical phases of the DFib and DYL models. 



B. Loss of the code property 



2. Ladder model 



We find that the non-Hermitian models are stable against 
local perturbations, and they satisfy a generalized code prop- 
erty. Keeping in mind that a non-Hermitian matrix has left 
and right eigenvectors, which in general are not identical, a 
local operator acts as a scalar multiple of an identity operator 
connecting the left and right ground state subspaces: 



(0 



(r) 



X(L)6 i: 



(11) 



Since only small linear dimensions are accessible to exact 
numerical diagonalization for the honeycomb lattice, we also 
consider a a quasi-one-dimensional ladder geometry consist- 
ing of rectangular plaquettes as shown in Fig. [5] The DFib 
and DYL models on this ladder geometry were introduced 
and solved in Refs.[T8landl4l respectively. Both models fea- 
ture topological phases with two (instead of four) degener- 
ate ground states, but are otherwise identical to the respective 
honeycomb lattice models. 
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FIG. 4. (color online) Scaling of the finite-size gap A(L) (in units 
of J p ) with linear system size for the Hermitian projector model 
jjhcrm Qn tWQ different lattice geometries: the honeycomb lattice 
with LxW plaquettes (top panel) and 2-leg ladder systems of length 
L (bottom panel). 



J r = sin 9 and J p = cos 9 , 

where 9 = corresponds to the unperturbed Hamiltonian. 
The phase diagrams as a function of 9 have been mapped out 
for both the DFib model 18 and the DYL model, 4 respectively. 

Directly probing the topological order in the DYL model 
and its Hermitian counterpart we show the lifting of their re- 
spective ground-state degeneracies in Figs. [6] and [7] when in- 
cluding a string tension. We find a striking qualitative dif- 
ference between these two models: For the DYL model the 
lifting of the ground-state degeneracy is exponentially sup- 
pressed with increasing system size - characteristic of a topo- 
logical phase. For the Hermitian model, on the other hand, we 
find a splitting of the ground-state degeneracy proportional to 
J r L. The linear increase with both system size and coupling 
can be easily understood by the different matrix elements of 
the string tension term on a single rung for the two degener- 
ate ground-states of the unperturbed model. Plotting the low- 
energy spectrum in Fig.[7]clearly shows that the two-fold de- 
generacy of the unperturbed Hermitian model arises from a 
(fine-tuned) level crossing. Similar behavior is found in the 
honeycomb lattice model (not shown). 

Considering the model in a wider range of couplings, as 
shown in Fig. [8] further striking differences between the non- 
Hermitian DYL model and its Hermitian counterpart are re- 
vealed: The DYL model exhibits two extended topological 
phases around 9 = and 8 = tt/2 (with two and four de- 
generate ground states, respectively), which are separated by 
a conformal critical point at precisely 9 C = 7r/4 as discussed 
extensively in Refs.|4]and[T8l In contrast, the Hermitian model 
jyherm gxjjjbits no topological phase anywhere, and the inter- 
mediate coupling 9 = 7r/4 does not stand out. 



die 
FIG. 5. Edge labeling for a plaquette of the ladder lattice. 



The quasi-one dimensional geometry allows to numerically 
diagonalize systems up to linear system size L = 13. The 
finite-size gap of the Hermitian model H herm is again found 
to vanish in the thermodynamic limit, showing a linear de- 
pendence on the inverse system size as shown in Fig.|4j3). To 
further demonstrate the fragility of these gapless ground states 
against local perturbations we add a string tension 18 



H pelt = J r 



rungs r 



(13) 



favoring the trivial label l(r) = 1 on each rung of the ladder. 
We parameterize the couplings of the competing plaquette and 



non-Hermitian DYL model 
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FIG. 6. (color online) Ground-state degeneracy splitting of the non- 
Hermitian doubled Yang-Lee model when perturbed by a string ten- 
sion (9 0). 
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FIG. 7. (color online) Ground-state degeneracy splitting of the Her- 
mitian model H hclul , the counterpart to the DYL model, when per- 
turbed by a string tension (6 7^ 0) (top panel). The slope of the 
splitting around the unperturbed model [0 = 0) is given in the inset 
(top panel) for different system sizes L. The bottom panel shows 
the low-energy spectrum, which clearly shows that the degeneracy at 
9 = is due to a level-crossing. 



IV. ABSENCE OF NON-UNITARY TOPOLOGICAL 
PHASES IN UNITARY MODELS 

So far, we have considered a specific set of Hermitian 
models constructed to have the same ground states as a non- 
Hermitian parent model and found that they no longer exhibit 
a topological phase. This raises the question whether this ob- 
servation points to a deeper principle, which we investigate in 
this section in rigorous mathematical terms. 

A. Galois Conjugate 

Let us now lay out the mathematical foundations as clearly 
as possible. The double DFib is isomorphic to a copy of Fib 
and its time reversal, DFib = Fib (g) Fib. Thus to Galois con- 
jugate DFib it is sufficient to define Fib e , then DFib e = 
Fib e ® Fib? . 

A theory such as Fibonacci can be defined using a set 
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FIG. 8. (color online) The low-energy spectra of the doubled Yang- 
Lee model (top) its Hermitian counterpart (bottom) for a wide range 
of coupling parameters. Data shown is for a ladder of length L — 8. 



of 6j-symbols braiding eigenvalues {-R^ } (not al- 

ways necessary), and some pivotal coefficients {e^ = ±1}, 
where k, I, m, n, a, b, c are anyon types (see chapter 4 of 
Ref. n~9T >. Because of gauge choices, there are many differ- 
ent sets for the same theory. If we fix a set of data, then 
we can define a number field K for a theory as the number 
field obtained from adjoining all numbers {i^ n } and {Ra C } 
to the rational numbers Q ( {e^ = ±1} are already in Q). 
The automorphisms of the number field K fixing Q form the 
Galois group of K, denoted as Gk- If 9 is an element of 
Gk, then by applying g to all data, we get a potentially new 
theory. We will call the new theory a Galois conjugate or a 
Galois twist. For the Fibonacci theory, the minimal number 
fields required for the Galois conjugation for both the alge- 
braic normalization and unitary normalization are worked out 
in Ref. 20 and needed below for the discussion of the projec- 
tors for code subspace property. For the algebraic normaliza- 
tion, the number field is the cyclotomic number field Q(^2o)> 
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where £jv = e 27ri / N , while for the unitary normalization, the 
number field is Q.(y/4>, £20)- 

It is known in general that a theory from quantum groups 
such as Fibonacci can always be defined within a cyclo- 
tomic field Q(£,n) for some N. For the Jones representa- 
tion with the algebraic normalization, this is done explicitly 
by Kuperberg 21 . To explain this, we digress briefly to some 
basic quantum topology. 

The Jones representation (and polynomial) may be con- 
structed from the Kauffman bracket: 



spaces via the automorphism of £{7(2)3: spinO = 

spin 3/2, spin 1/2 = spin 1, spin 1 = spin 1/2, spin 3/2 = 
spinO. Similarly, X(n ■ 2, 0, e 4 "/ 5 ) X(n • 1, 0, e 47 ™/ 5 ) as 
Hilbert spaces of mixed sign. The braid group actions (Galois 
conjugates of the Jones representation) are, of course, also 
identical. For t = e 4 * 1 / 5 , the loop value is d = — e 27 ™/ 5 — 
e -2W5 = (f) = ii^S, the golden ratio. 

We use trivalent graphs with the Kauffman vertex normal- 
ization: 




t indeterminate. 

The skein space W(n ■ 1, 0, t) is the vector space of formal 
linear combinations of arc matchings (i.e. skeins) of an even 
number, n, of fixed points on the top of a square (the arcs are 
imbedded in the interior of the square and "0" means no points 
marked on the bottom of the square). An n-braid b acts on W, 
by gluing b onto the top of the square and resolving crossings 
by the above rule. The trick is to define each braid genera- 
tor as t 1 / 4 times the geometric crossing. Additionally, each 
skein is "even" or "odd" according to whether a 2-coloring of 
the skein complement (starting with the bottom of the square 
being colored white) has an even or odd number of black re- 
gions. All black skeins should be multiplied by a factor of 
t 1 / 2 . This results in a basis for W and a rescaling of the ac- 
tion so that this 'W-representation" is defined over the field 
Q[i\. The famous quantum representations of Jones at level k 

2tti 

will be quotients of the W^-representation for t = e k + 2 . Note 
that the rescaling cannot affect the density of the projective 
Jones representation, which will be important shortly. 

In the case at hand, Fib, t = e 2Trl ^ 5 , and the Galois conju- 
gate theory Fib e is obtained by replacing t by t' = e 47 ™/ 5 . 

The skein space W(n •]_, 0, t) carries a natural bilinear form 
( , ) obtained by doubling the square (thought of as a disk) 
along its boundary and evaluating the union of the two skeins 
as a scalar using the above Kauffman relations. When \t\ = 1, 
the form (A, B) is Hermitian. If further Hs a root of unity, 
then this form has a singular subspace. X(n ■ 1,0, t) is, by 
definition, the finite dimensional Hilbert space obtained by 
annihilating this kernel. The Hermitian form {A, B) is non- 
singular on X and the braid group B n acts. 

When t — e k + 2 , this quotient action is the Jones represen- 
tation associated to SU(2)k, whose trace leads to the Jones 
polynomial evaluated at t. For t = e 27Tl /(. k + 2 )^ the Hermi- 
tian form (A, B) is positive definite. For other roots of unity, 
(A,B) may be of mixed signs (p,q), p 7^ 0, q 7^ 0. This 
happens in particular for t' — e 47T W 5 when n > 4 as we now 
check. 

Well-established conventions in mathematics and physics 
lead to two different ways to label the particle types in 
SU(2)k'. one by the spins of the irreps, and the other by the 
dimensions of irreps minus one. Unless we speak explicitly of 
a spin label, as in the next paragraph, the labels of particles in 
this section are by the dimensions minus one, which are twice 
of the physical spins. 

Note that X(n-2, 0, e 27 ™/ 5 ) ^ X(n ■ 1, 0, e 27 ™/ 5 ) as Hilbert 



1 




where TT is the Jones-Wenzl projector 



and 1/2 and 1 are the spins of the quasi-particles, to write an 
orthogonal basis for X(A ■ 2, 0, t) (The ".2" after n indicates 
the second, i.e. spin 1, nontrivial particle type). 



WW 




WW 










HO 



is a positive number. Now consider an orthogonal basis ele- 
ment 



We have 




Its self pairing is negative: 




a negative number where d is the loop value W= and the 6- 
symbol = — <f> (simplified from the formula in Figure 19 in 
Ref. |22| |. It follows that for n > 4 the Hermitian structure 
on X(n ■ 2, 0, e 47 ™/5) ^ X {n ■ 1, 0, e 47 "/ 5 ) has mixed signs. 
The corresponding braid representations for DFib e also have 
mixed signs when n > 4. 

The doubled Fibonacci theory DFib, as with all topological 

phases, has the code property: the composition G ^> H —> 

p 

H — > G is multiplication by some scalar Al € C whenever 
L is a (sufficiently) local operator (P = inc^ is the Hermitian 
orthogonal projection to the Levin- Wen ground state G, see 
Ref. |8). In general P is Hermitian, but for Fib P is actually 
real symmetric. 



8 



In the Levin-Wen Hamiltonian scheme, there are two kinds 
of terms: the vertex type A v for each vertex v, and a plaque- 
tte type B p for each plaquette p. The code space G is the 
common eigenspaces of all local operators {A v } and {Bp}. 
With the standard choices of basis, the vertex terms A v are 
matrices with entries O's and l's, while the plaquette terms 
B p are matrices with entries given by products of 6j-symbols. 
The algebraic constraints defining the code subspace trans- 
form under a Galois conjugation. Hence solving the Galois 
conjugates of the constraints defining G, we obtain G e and 
replace the matrix P with its Galois conjugate P s . Clearly 

G Q ^> H A H ^ G Q is multiplication by (A i(s -i,) e 

since L is local if and only if i/ e ) is local. Thus G e re- 
tains a "code" property but with respect to a non-Hermitian 
projector P G . In the symmetric normalization of Yang-Lee, 
P s is complex symmetric with eigenvalues and 1 (since 
(P e ) 2 = P G , however (P e )t = P<3 f P Q ). 

We close this paragraph by noting that our proof uses the 
mixed signatures in the Jones braid group representations as 
shown above, and as such applies only to the algebraic nor- 
malization. For the symmetric normalization of the DYL the- 
ory, the Jones representation spaces have either positive or 
negative definite inner products, but the mixed sign in the al- 
gebraic normalization will be sufficient to prove our theorem 
for any choice of normalization. 



B. Lieb-Robinson bounds and local unitary evolution of a 
ground state under changes in the Hamiltonian 

Our proof of absence below will be based on a contradic- 
tion of the above result for the Galois conjugated theory with 
local unitary evolution of a ground state under local changes 
in a Hermitian Hamiltonian. This local unitary evolution can 
be proven for all Hermitian Hamiltonians that satisfy Lieb- 
Robinson bounds. 

Lieb-Robinson bounds are a mathematical way of express- 
ing the physical fact that in local lattice Hamiltonians there 
is some upper bound to the velocity of excitations. These 
bounds can be proven for a wide-range of Hamiltonians, in- 
cluding what would be colloquiually referred to as "Hamil- 
tonians with finite-range interactions" or "Hamiltonians with 
exponentially decaying interactions". 

For a precise statement of conditions under which Lieb- 
Robinson bounds can be proven, we follow Ref. |23]where one 
sufficient condition is given as follows (see also Ref. l24t . We 
consider lattice Hamiltonians, and use ... to label sites of 
the lattice, with some metric dist(z, j) on the lattice. We use 
X, Y, Z, ... to label sets of sites of the lattice. Let the Hamil- 
tonian H be written as 



(14) 



where the operators Hz are supported on sets Z (an operator 
is said to be supported on a set Z if it can be written as a 
tensor product of an operator on the degrees of freedom on 
set Z with an identity operator on the remaining degrees of 



freedom). Assume that the following condition holds for all 
sites i, 

\\H x \\\X\exp[fji diam(X)] < s < oo, (15) 

for some positive constants /i, s, where diam(X) denotes the 
diameter of set X, and \X\ denotes the cardinality of X, and 
denotes the operator norm. 
Then, 23 Eq. (15i implies the following Lieb-Robinson 
bound for Hermitian Hamiltonians. For any operator O, we 
use 0(t) to denote the Heisenberg time evolution of the op- 
erator: 0(t) = cxp[iHt}0 exp[— iHt]. Let Ax , By be op- 
erators supported on sets X, Y, respectively. Then there is a 
constant vlr depending only on s, /i such that for t real with 
\t\ < dist(X, Y)/vlr, we have 



[A x (t),B Y ]\\ < 



vlr\A 



g{l)\X\\\A x \\\\B Y \ 



(16) 



where / = dist(X, Y) and g(l) decays exponentially in /. 
Given that Ax(t) has small commutator with all operators 
By with sufficiently large distance dist(Jf, Y), this implies^ 
that the operator Ax (t) can be approximated by an operator 
A l x (t) which is supported on the set of sites within distance 
I — VLn\t\ of the set X up to an error in operator norm which 
is bounded by ^f^-g{l)\X\\\A x \\. 

Definition IV.l. We say that a Hamiltonian is a Lieb- 
Robinson Hamiltonian (or that it obeys a Lieb-Robinson 
bound) if a bound of the form Eq. ( [To} holds for some v^r and 
some exponentially decaying g(lj. A parameter dependent 
family of Hamiltonians H s uniformly obeys a Lieb-Robinson 
bound if for some vlr and g(l) the bound Eq. ( 16 1 holds for 
all s. Such a family is called "uniformly LR". 

We also want to define what it means for a Hamiltonian to 
have multiple ground states and a spectral gap. Note that it 
is common practice in physics to refer to a system, such as 
a fractional quantum Hall system which has three low-lying 
states with an exponentially small splitting between them and 
then a gap to the rest of the spectrum as having a "degenerate 
ground state", even though the non-vanishing splitting means 
that the lowest eigenvalue is in fact non-degenerate. Our def- 
inition will reflect this usage, as we will not require that the 
states that we refer to as "ground states" be degenerate. All 
we will require is that the "ground states" be separated from 
the rest of the spectrum by a gap. 

Definition IV.2. A Hamiltonian has n ground states and a 
spectral gap AE, if E n -i + AE < E n where the eigenvalues 



of the Hamiltonian are E 0} E 



with E a < E 1 < 



family of Hamiltonians H s has n ground states and a uniform 
spectral gap AE if £J n _i(s) + AE < E n (s) for all s, where 
the eigenvalues of H s are Eq(s) < -E'i(s) < ... 

Note that we did not require in the above definition that the 
splitting E n -\ — Eq between the different ground states be 
small. For all the systems we are concerned with, this split- 
ting will turn out to be small, but since it is not required to 
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be small for lemma |IV.3| we do not include this in our def- 
inition (in some applications of Lieb-Robinson bounds, the 
splitting between different ground states is important, but we 
don't need it here). 

The next lemma expresses how the ground states evolve un- 
der changes in the Hamiltonian. We consider some parameter 
dependent family of Hamiltonians, H s , for < s < 1, and 
imagine this family as describing some path from an initial 
Hamiltonian at s = Q to some final Hamiltonian at s = 1. 
Stated roughly, this lemma shows that if the Hamiltonian is 
gapped and local, then the change in the ground state under a 
local change in the Hamiltonian can be expressed by a local 
operator acting on the ground state. 

Lemma IV.3. Let H s be a uniformly Lieb-Robinson family of 
Hermitian Hamiltonians, for < s < 1, with H s differen- 
tiate with respect to s, such that d s H s is supported in a disk 
X of radius R and such that for some J, and \\d s H s \\ < J for 
all s. Let H s have uniform gap AE. Let P(s) denote the Her- 
mitian projector onto the ground state subspace of H s . Then, 
for any I, there exists a family of unitaries U s supported on 
the set of sites within distance I of X such that 

\\U 3 P(0)U^-P(s)\\ < const.^(e X p(-/AB/2w iR )+ 5 (Z) 

(17) 

Proof. The proof largely follows previous results on quasi- 
adiabatic continuation and is given in Appendix [A] for com- 
pleteness. □ 

We make a few remarks. First, note the appearance of g(l) 
in the lemma above. For a Hamiltonian with exponentially 
decaying interactions, g(l) will decay exponentially in I, but 
for a Hamiltonian with bounded range interactions, g(l) will 
decay faster than exponentially in I, and the error will be dom- 
inated by the term exp(— IAE /2vlr)- Further, in the case 
of exponentially decaying interactions, the length scale over 
which g(l) decays will be set by the decay scale of the inter- 
actions in the Hamiltonian, i.e., by the microscopic details of 
the interaction rather than the magnitude of the spectral gap. 

Note that in the lemma above, a bound on ||9 s Jf s || appears; 
that is, the bounds depend upon how rapidly the Hamiltonian 
changes along the path. To give a physical explanation of why 
this appears, consider dragging an anyon along some path. 
Suppose we move the anyon a distance L. Then, since we 
always scale the path length to unity (that is, s ranges from 
to 1), the "velocity" at which the anyon moves along the 
path is proportional to L. Thus, for larger L, we are moving 
the anyon more rapidly along the path, and so ||<9 s _ff s || will be 
larger; thus, in a sense the appearance of ||9 s -ff s || is really a 
way of measuring the distance we drag the anyons. Thus, it 
is worth restating the result in a re-scaled way: suppose that 
we drag an anyon a distance of order the disk radius R. Then, 
typically we will have ||9 s iJ s || < J for a J of order R. So, 
given that the error in Eq. ( 17 1 is exponentially small in I, in 



Also, note that if d s H s is approximately supported in X, in 
that it can be approximated, up to exponentially small error, 
by an operator supported in X, then we can derive a similar 



bound to Eq. ( 17 1 which will involve the error in approximat- 
ing d s H s . We omit this case. 

Finally, in the case that the Hamiltonian H s is a sum of 
commuting terms with bounded-range for all s, the Lieb- 
Robinson velocity is zero. In this case, it is possible to show 
that for sufficiently large I, the error \\U s P(0)U]. - P(s)\\ is 
exactly zero. 



C. Proof of Absence 

Examined in detail, the ground state manifold G e and the 
projector P Q which defines it depends on: 1) the number and 
location (V) within the 2-sphere S 2 of the anyons, 2) the anyon 
particle type — a kind of boundary condition, and 3) the (possi- 
bly non-unitary) trivalent vertex normalization / : L 3 — > C\0 
or gauge choice, L being the label set. For Fib, DFib, and 
their Galois conjugates and time reversals (represented by"), 
/ is always symmetric and satisfies a consistency relation with 
the F-symbols: suppose {Ff? } are new 6j symbols from 

{ F inm\ b y a S au § e chan g e {/(a, b, c)}, a,b,c£ L, then 



F, 



ijk pijfc fU,k,n)f(i,n,l) 



l:nm 



l ' nm f(i,j,m)f(m,k,lY 



Except that it would unpleasantly cluster the notation, we 
should write G n r ^ and P^ r ^. The detailed position T of 
the anyons within the lattice model is important to us since 
our proof will work with the entire "braid groupoid" B n . In 
fact, we treat T as a continuous variable on a compact space 
of 2n (real) dimensions. This moduli space of anyon po- 
sition is compact since distinct anyons are not permitted to 
closely approach. The elements of B n are oriented paths of 
n-distinct (marked and framed) points in M. 2 which compose 
only when end points match. B n represents in a large but fi- 
nite dimensional Hilbert space H of microscopic degrees of 
freedom on S 2 , the north pole serving as a standard oo for 
M 2 . The vertex normalization / is also important within the 
proof. As we have already seen, the symmetric normaliza- 
tion yields a TQFT with all definite Hilbert spaces (though 
some are positive-definite and others are negative-definite). 



such a case it suffices to choose I logarithmically large in R 
in order to make the error of order 1. For an I of order R, the 
error will be exponentially small in R. 



The proof of Theorem IV.5 requires as a "kernel" a single 
Hilbert space on which a non-singular form of mixed signs is 
preserved by B n . With this kernel in hand, the proof actually 
covers all vertex normalizations /. 

Definition IV.4. We call an operator C range r if it is sup- 
ported on a ball of diameter r. Also, we use the same term 
for sums of such operators. Similarly, an operator is called 
weakly range r (in either sense) if it is range r up to expo- 
nentially small corrections. We say that an operator is short 
range if it is supported on a ball of diameter small compared 
to system size. 

We say an operator O is a local normalizer iff there is some 
constant c which is small compared to system size such that 
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OLO~ x is range r + c whenever £ is range r. We say that an 
operator is a weakly local normalizer iff there is some constant 
c which is small compared to system size such that OCO^ 1 is 
weakly range r + c whenever C is range r. 

A uniform family of (weakly) local normalizers 0\ is a 
parameter dependent family of operators such that O^CO^ 1 
is (weakly) range r + c whenever C is (weakly) range r, with a 
uniform bound on the exponentially small corrections and on 
the constant c, and such that whenever |A — A'| < 0(1), the 
product 0\0~^} is a product of at most 0(1) operators which 
are all (weakly) range r and are not necessarily the same, for 
some r which is 0(1). An example of a local normalizer is 
a finite depth quantum circuit of invertible (not necessarily 
unitary) local operators. An example of a uniform family of 
local normalizers is a family of finite depth quantum circuits 
of invertible local operators, such that an 0(1) change in the 
parameter changes only 0(1) different operators in the circuit; 
for the applications we have in mind, one should imagine that 
the parameter A refers to different anyon positions and that 
changing A changes the circuit only near the anyon positions. 

In the definition of weakly local normalizer, it will be im- 
portant to define how we quantify the error term in the approx- 
imation by a bounded range operator. The natural way to do 
this would be to require that the error term be small in operator 
norm compared to the operator norm of OCO^ 1 . However, 
for technical reasons, for use later we will be interested in 
a what we call a g.s. -weakly local normalizer (g.s. stands for 
ground state). In this case, we consider certain operators M(i) 
which have the property that M(i) is bounded range and ex- 
actly maps the ground state subspace of some non-Hermitian 
Hamiltonian to the ground state subspace of some other non- 
Hermitian Hamiltonian, with M(i)^ M (i) exactly preserving 
the ground state subspace of the first non-Hermitian Hamilto- 
nian and having its ground state expectation value equal to its 
norm. Then, we require that the error term be small in oper- 
ator norm compared to the norm \0 A(i+i)iW for ip 
in the ground state of some other Hermitian Hamiltonian (this 
ground state subspace is obtained by applying O to the ground 
state subspace of the Hermitian Hamiltonian). Note that if O 
were an isometry, then the norm \0\^ + i)M(i)0^^tp\ would 
equal the norm of M and so this would reduce to the more 
natural definition. Note also that any local normalizer is a 
g.s. -weakly local normalizer. 

Theorem IV.5. Fixing the number n > 5 and particle type 
t ® t of D Fib anyons on S 2 and any vertex normalization f 
there can be no continuous uniform T-family of (g.s. -weakly) 
local normalizer operators Op-' H — > W, so that OrG^ r ^ is, 
for all anyon positions T the ground state manifold of a uni- 
formly Lieb-Robinson and uniformly gapped family of Hermi- 
tian Hamiltonians H(T) defining a topological phase [see Eq. 

07- 

Proof. The theorem uses the notation of reference 21 to 
describe the anyons in DFib. For now, fix the algebraic vertex 
normalization A = /. Below we may suppress T and / from 
the notation when they play no role. 

Suppose 0p exists, then0pGp is a family of code sub- 
spaces and for T near T' the subspaces are connected up to ex- 



ponentially small disc repan cy by a local unitary £/r,r' (these 

Fib? 



are the U s of Lemma 



IV.3 



Fib / 



Writing DFiby = nu^ « nuj 
(one may think DFib describes a bilayer), let us recall a the- 
orem stated in|2T|for the right hand factor Fib^, where / is 
the algebraic normalization (Note: while DFib^ is a theory 
of string -nets on the surface S 2 , with boundary conditions at 
anyons, Fib^ is the corresponding string -net theory^ in the 3- 
ball with boundary S 2 . Thus the function / gauging vertices 
acts compatibly in both theories.) 

Now, according to Ref.BTl corollaries 1.2.4 and 1.2.6, for 
n > 5, the Jones representation p on the topologically defined 
Hilbert space Vr of ground states for Fib^ is (analytically) 
dense in a noncompact special unitary group (preserving a 
Hermitian metric of mixed signs) SU(p,q) := SU((X,n ■ 
1, 0, e 4 ^)) £ SU((X, n ■ 2,0, e 4 ™/ 5 )). 



Recall that in subsection IV A we confirmed that the loop 



and that with the 



value for a closed string in DFib^ was 

algebraic vertex normalization the signs of the B n -invariant 
ground state "Hilbert" spaces are indeed mixed: p > and 
q > 0. 

We need to formulate a lemma regarding the following con- 
cept. 

Definition IV.6. An ambient groupoid representation is a 
functor from a groupoid to the category of subspaces of a fixed 
Hilbert space, and linear transformations on H carrying one 
image subspace to another. 

Thus to objects a, b of the groupoid we assign spaces A C 
H, B C H and to a morphism a — > b a unitary map H — > 
H carrying A to B. In our case, the groupoid is B n with T 
being the objects and motions of anyons being the morphisms. 
The subspaces are the respective ground states for H^ w alg — 
the Galois conjugated Levin- Wen Hamiltonian with algebraic 
vertex normalization. 

Lemma IV.7. Let {A,}, i € T be the set of code subspaces of 

fixed Hilbert space H and B = {T, {— >}} a groupoid. Sup- 

q 

pose for a generating set of morphisms in B, a — » b, there 
are invertible local operators Cq : H — » H with LA = B 
(A=image a and B=image b). Then this data determines a 
unique projective ambient groupoid representation . That is, 
all compositions commute with the functor up to multiplica- 
tion by a nonzero scalar. 



C e \ A 



C<j,\i 



Proof. Consider A — ^ B C, C e and £ : H -> H 
being the local operators carrying A to B and B to C respec- 
tively. The composition ■ Ce is local and ■ £e(A) = C. 
By the code property there is a scalar A so that Tl^C^-Cg)" 1 ■ 
(Ctf, ■ Ce) ■ incA = A • Id, where A ^ since all morphisms C 
in the representation are invertible. (n^ denotes the appropri- 
ate projection to A corresponding to its code property.) Thus 
(projectively) the homomorphic property, restricted to the sub- 
spaces {Ai}, i 6 r, is redundant when these subspaces have 
the code property. 

For uniqueness consider two possible representations C(9) 
and C'(6). {C')- 1 C : H -> H is local and carrying A to 
itself. Thus n j 4(£g) _1 £g • inc^ : A — > A is multiplication by 
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some scalar A which, again by the invertibility of Cg and C' g , 
is non-zero. □ 

By assumption the collection {0(Gp)} is a code with re- 
spect to the usual Hermitian projection. Using the U s of 
Lemma [IV. 3 1 as the generating set of morphisms, Lemma [IV.7| 
builds a projective representation p : B„ x H H of DFib 
up to exponentially small errors which we can neglect for the 
moment but return to shortly. We may think of this represen- 
tation as the result of (quasi-)adiabatic evolution of Gp in- 
side the finite dimensional microscopic Hilbert space H. In- 
tuitively, braiding might be realized by building and slowly 
moving a potential trap term added to the Hamiltonian. For 
plaquette excitations (e.g. t (g) r) such a trap could have the 
rough form H trdp — arctan(St)B p +(ir—arctan(8t)B p >+e(j* 
for plaquettes p and p' separated by edge. Such family of 
Hamiltonians will adiabatically braid the anyons. Formally, 
however, we have posited, for contradiction, the family H (T) 
and this is all we need. 

provides local unitaries intertwining Gp and 



Lemma 



IV.3 



G^ 2 , and so p is obviously a (projective) unitary representa- 
tion with respect to the standard positive definite Hermitian 
structure on the space H of microscopic degrees of freedom. 
Thus p preserves ordinary lengths and angles (as measured in 
H) and the complex structure (multiplication by i) of H as 
well. So this p manufactured from Lemma IV.7 looks geo- 
metrically quite distinct from a second representation of B n , 
p' := 0(p <g> p*)0~ 1 . As explained above (also see refP^I), p 
and p* also act on the same collection of string-net spaces as 
p and the uniqueness clause of Lemma |IV.7| yields a projective 

proj 

isomorphisms p = p' . 

We learned from Ref. [2TI that p is dense in SU(p,q),p > 
0, q > 0, and that X(n ■ 2, 0, e^s 1 ) is the fundamental repre- 
sentation uj\ of SU (p, q). oj\ preserves a form of mixed signs 
and thus will distort not only Euclidean length but Euclidean 
angles as well by an unbounded amount (For example con- 
sider the effect of boosts in 0(1, 1) C 17(1, 1) on Euclidean 
angle). This difference, that p and p' preserve forms of differ- 
ent signatures, excludes the existence of a (g.s. -weakly) local 
normalizer O transforming {G^} to the ground state spaces 
ofP(A). 

We now consider in more detail the exponentially small er- 
rors that we have neglected. First some preliminaries. We will 
need the following lemma which provides a useful corollary 
of the disk axiom: 

Lemma IV.8. Let X be some set and P some projector such 
that for any operator A supported on X there is a scalar z 
such that 



\PAP- zP\\ < e\\A\\ 



(18) 



for some sufficiently small e. (In typical applications, we have 
in mind that P is the projector onto the ground state subspace 
of some system, X is some set of small diameter, and the above 
equation encodes a soft form of the disk axiom for that theory). 
Then, for any operator O supported on set X there is a scalar 



w such that 

\\POP - wP\\ < GVejKl - P)OP\\ \\PO(l -P)\\, (19) 

for some constant C of order unity ( the constant C is indepen- 
dent of e, so long as e is sufficiently small). 

Proof. First, assume that O is Hermitian. Consider the opera- 
tor U = exp(UO), fort real. Since \\U\\ = 1, \\PUP-zP\\ < 
e for some z. We can expand PUP in a power series in 
(1 - P)OP giving 

PUP = Pexp(itPOP)P + 0(t 2 \\PO(l - P)|| 2 ). (20) 

Suppose POP = w + A, for some traceless operator A and 
scalar w. Pick t = 2e/||A||. Then, 

PUP = Pw+2ieP^P+e>(e 2 +e 2 ||PO(l-P)|| 2 /||A|| 2 ). 

(21) 

So, the distance from PUP to the closest scalar multiple of 
P is at least 2e - (D(e 2 + e 2 \\PO(l - P)\\ 2 /\\ A|| 2 ). If || A|| 
is sufficiently large compared to y/e\\PO(l — P)\\, then the 
term 0{...) is small compared to the leading term, assuming e 
is sufficiently small (we need e sufficiently small compared to 
unity so that term 0(e 2 ) is small). However, this contradicts 
the assumption that PUP is within distance e of some scalar 
multiple of P. 

Now, consider the general case that O need not be Hermi- 
tian. Add an additional spin- 1/2 degree of freedom to the 
system and consider the Hermitian operator O = (O — w) <S> 
(7+ + (O* — w) ® <t~ , where a + ,a~ are the raising and lower- 
ing operators for that spin and w is chosen so that P(0 — w)P 
is traceless. The assumption ( 18 1 for the original theory im- 
plies that for the system with the additional spin- 1/2, for any 
operator A acting on set X and on the additional spin- 1/2 that 



\PAP-P®Q\\ <4e||A||, 



(22) 



for some 2-by-2 matrix Q acting on the additional spin (to 
show the above equation, expand PAP as a sum of four prod- 
uct operators, one operator in the product acting on X and the 
other on the added spin, and apply Eq. ( [T8| to each term in the 
product). Construct U = exp(iiO). Perturbatively expand 
the a + component of PU P for the same t as before. This is 

2 i eP^P®(7 + +C'(e 3 +e 2 || J PO(l-P)||||(l-P)OP||/||A|| 2 ) 

(23) 

We again get a contradiction as in the Hermitian case. □ 

The representation p gives a mapping from braids to ma- 
trices. We will use M to refer to such a matrix. We can also 



construct a matrix U by taking the U s of lemma IV.3 for the 
corresponding braid and projecting into the grounds state sub- 
space. So long as the length of the braid is smaller than some 
quantity growing exponentially with the linear size of the sys- 
tem, the matrix U will be approximately unitary (the error 
arises from leakage out of the ground state subspace in lemma 
|IV3] We claim that 
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Lemma IV.9. 



So, 



\U - zGMQ- 



« 1 



(24) 



for such braids for some scalar z. Indeed, the difference in 
norms is exponentially small in system size. 



Proof. This difference follows from the disk axiom: we de- 
compose the braid into n segments; taking n of order L, each 
segment moves the anyons only a short distance. In the non- 
unitary theory, we can construct short range operators which 
move the ground state subspace before the i-segment to the 
ground state subspace after that segment. Call these operators 
M(i), so that M is equal to the projection of M(n)...M(\) 
into the ground state subspace. This operator M (n)...M(l) 
exactly preserves the ground state subspace of the Galois con- 
jugated theory and OMO^ 1 exactly preserves the ground 
state subspace of the unitary theory. We use P(i) to denote 
the projector onto the ground state subspace after the i-th seg- 
ment, and we use 



M'{i) = O m M{i)Oll 



(25) 



where Oa(i) is the operator 0\ corresponding to the position 
of the anyons after the i-th segment, so that M(i)P(i — 1) = 
P(i)M(i) and P(n) = P(0). In the unitary theory, let U(i) 
denote the unitary matrices from lemma |IV.3| for the motion 
along the i-th segment, so that U is equal to the projection 
into the ground state subspace of U (n)...U(2)U(l). Note that 
U(n)...U(l) preserves the ground state subspace up to ex- 
ponentially small exponentially small error and the matrices 
U(i) are bounded range. Thus, 

\\U - zOMO~ x \ (26) 
= \\P(0)U(n)U(n - 1)...C7(1)P(0) 

-zP(0)M f (n)M'(n - 1)...M'(1)P(0)|| 
« \\P(0)U(n)P(n - l)U(n - l)P(n - 2)...U(1)P(0) 

-zP(0)M'(n)P(n - l)M'(n - l)P(n - 2)...M'(1)P(0)||, 

where the left side of this approximate equality differs from 
the right by inserting additional factors of P(i — 1) after every 
occurrence of U(i) or M'(i) on the right-hand side. We have 
P(0)M'(n)M'(n - 1)...M'(1)P(0) = P(0)M'(n)P(n - 
l)M'(n - l)P(n - 2)...M'(1)P(0), but the error in Eq. ((26 1 
occurs because U(i)P(i — 1) is only approximately equal to 
P(i)U(i — 1) (the difference is exponentially small in system 
size). Note that on the right-hand side of this equation the 
matrices act in the ground state Hilbert space, while on the 
left-hand side they act in the full Hilbert space. To bound the 
right-hand side of Eq. (|26|, it suffices to show that, for all i, 



\\P(i)U(i)P(i - 1) - z(i)P(i)M'(i)P(i - 1)|| (27) 
is exponentially small for some scalar z(t). To see this, set 



\\P(0)U(n)P(n - l)U(n - 1)...P(1)E7(1)P(0) - (29) 
P(0)M f (n)P(n - l)M'{n - 1)...P(1)M'(1)P(0)|| 

< \\P(0)U(n)P{n - 1)...P(1)C/(1)P(0) - 
P(0)U(n)P(n - l)z(n - l)M'(n - l)...z(l)M'(l)P(0)|| 
+ \\P(n)U(n)P(n - 1) - P(n)z(n)M'(n)P(n - 1)|| 

xY[\\p®M'®p(i-m 

i=l 

< ... 



n^) 



(28) 



Given this norm estimate (27 i, then since P(i)U(i)P(i — 1) 
is an approximate isometry from the range of P(i — 1) to 
the range of P(i), the matrix P(i)z(i)M'(i)P(i - 1) is 
also such an approximate isometry, so the product of norms 
\\P(i)M'(i)P(i - 1)|| above is bounded. 

So, we must bound Eq. |27] >. Since Ui) is an approx- 
imate isometry, it suffices to bound \\P(i — 1) — P(i — 
l)U(iy ' z(i)M'(i)P(i - 1)||. At first sight, this seems to fol- 
low immediately from the disk axiom: since U (iY M' (i) is 
short range, or at least approximately short range (note that 
for M'(i) this follows by the definition of a family of lo- 
cal normalizers, but see the next paragraph for a more care- 
ful treatment of error terms), by the disk axiom it is close 
to a scalar when projected into the ground state subspace. 
Hence, choosing z(i) to be the inverse of this scalar, the de- 
sired result seems to follow. However, there is a complica- 
tion: suppose P(i — l)U(iyM'(i)P(i — 1) is within some 
distance e of z(i)~ 1 P(i — 1) for some z(i); then, we bound 
\\P(i - 1) - P(i - l)U{i)U(i)M'(i)P(i - 1) < e\z(i)\. 
Hence, if z(i) is large, the resulting error can be large even if 
e is small. This is why we will need the lemma ( |IV.8| l above. 

By definition of g.s. -weakly local normalizer, the operators 
M'{i) can be approximated by operators that are short range, 
up to an error that is small compared to \M'(i)ip\ for all ip in 
the ground state subspace P(i — 1) with \tp\ — 1. Since U(i) is 
approximately unitary and an approximate isometry between 
two ground state subspaces, this means that U(iyM'(i) can 
be approximated by an operator 0(i) that is short range, up 
to an error that is small compared to \U(iyM'(i)ip\ for ip in 
P(i — 1). (Note that if O is a local normalizer, then M'(i) 
already is short range so we can take 0(i) = U(iyM'(i) in 
that case.) So, for the 0(i), ||(1 - P(i - l))0(i)P[i - 1]| is 
small compared to |0(i)V>| for all ip. Applying lemma ( |IV.8 1, 
this means that P(i — l)0(i)P(i — 1) is close to z(i)P(i — 1), 
for some z(i) up to an error that is small compared to z(i). □ 

We can find a braid such that the corresponding matrix M 
is diagonalizable, and with the ratio between its largest and 
smallest eigenvalue being at least 2 in absolute value^. This 
means that the ratio between the largest and smallest eigen- 
value of zOMO^ 1 is at least 2 in absolute value. However, 
zOMO~ x is close to a unitary matrix. All eigenvalues of a 
unitary matrix are on the unit circle in the complex plane. Fur- 
ther, a small perturbation of a unitary matrix leaves all of its 
eigenvalues close to the unit circle, so zOMO^ 1 must have 
all of its eigenvalues close to the unit circle contradicting the 
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assumption on the ratio of eigenvalues. 



But 



Now we remove the condition on the vertex normalization 
or "choice of gauge". As noted as above the gauge choice is 
a function / : L 3 — > C\0. Although this is not crucial, it is 
pleasant that in the Fib case / is identically 1 except for tak- 
ing value /(t, t, t) = A on the essential trivalent vertex. Re- 
calling the A v and B p terms in the Levin-Wen model Hamil- 
tonian, the fusion rule terms, A v , commute with and do not 
depend on /. The effect of / on B p may be computed (using 
the compatibility of gauge choice and the F-matrices used to 
construct B p ): 



inc^s oFloO lj/ o£oO'oF\omc n e = inc n goO ,ji oCoO'omc, 

(35) 

So we find that for the change of variables O', C acts on 
O' o inc G s . Varying O over all local normalizer operators pro- 
duces a local normalizer O' = O o F x . Thus for every pos- 
sible local normalizer operator O 1 there is a local C : H — > H 
acting on O'G x . 



This completes the proof of Theorem |IV.5| by removing the 
hypothesis of algebraic vertex normalization. □ 



B P: x = F\ ■ B p>a i g , ■ F x , 



(30) 



where Fa is the "relative fugacity matrix", Fx : H — > H . 
Fx is diagonal in the string-net basis {k} and the entry Fk,k 
is simply \#( k \ where #(fc) denotes the number of (r, r, r)- 
vertices in the k th string-net. 

As a concrete example of the above gauge dependence of 
B„ in the Fib case 



pTTT 

T;sj/m 



1 



\-2) F T-alg l A 2 



1 



(3D 



By inspecting the two F-matrices in the two normalizations, 
we have A 2 = i0 5 / 2 . Given a trivalent graph 7, and a labeling 
of its edges fc 7 . Let {& a lg.,k^} and {^ sym ..k-,} be the two 
basis of the Hilbert space (g) ee7 C* ( ^ L - ) of labeled graphs. Then 



sym.,K^ 



A #( ^)* Q / fl .,fe 7 - Suppose 



Bp^alg.^alg..k,, = / „ -Bp,atg.,fe;,k r ^ alg.,k^, (32) 



noting that gauge change and recoupling are commutative, we 
obtain 



B 



(33) 

Thus Br.sym has the claimed conjugated form. Similarly for 



any vertex fugacity A, Eq. (30 1 holds. Thus for a general 



relative fugacity A, Q r ,X = F x Gr,ai g . 

Observe that while Fx is not local (in the sense of hav- 
ing supported on a disk of bounded radius), both Fa and 
Fr 1 are implemented by a depth 1 invertible circuit, Fx — 
n it (AII rTT + (1 — II rrT )). Thus the general local normal- 
izer operator O may be written as O = O' o Fx, where O' 
is also local normalizer; given X,0, and O' determine each 
other uniquely, we have just shown that for all local normal- 
izer O there is a local L so that C acts on OG^ lg , i.e. that 



o 0t o C o O o inc G s ^ scalar. 



(34) 



It follows that for all O' there is an C (identical to C above) 
so that: inc G s o F\ o O't o C o O' o Fx ° inc G s ^ scalar. 



Theorem|IV.5 immediately implies the following corollary 



Corollary IV.10. Let {Gp ^} be the ground state manifolds 
for the Galois conjugated Levin-Wen Hamiltonian Hf, ^ for 
n > 5 t (E> T-anyons on the 2-sphere S 2 with positions T 
and any vertex normalization f, within a larger Hilbert space 
H of microscopic (lattice) degrees of freedom. There can be 
no continuous uniform family of local normalizer operators 
Or so that {OrHf, ^Cp 1 } are uniformly gapped uniformly 
Lieb-Robinson Hamiltonians determining topological ground 
states {OrGp A, in the sense ofTQO-P 2 (i.e. satisfying the 
code property (f7) ). 

Although we have concentrated the discussion on the Fib 
TQFT, its quantum double, and their Galois conjugates, the 
proof requires only two ingredients: 1) finding pairs of Galois 
conjugate theories (with choice of vertex gauge) one of which 
is unitary (for the Hilbert space of a sphere or plane with fixed 
anyon content) and one of which is unitary with respect to a 
Hermitian metric of mixed signs (p, q),p > 0,q > 0, and 2) 
establishing denseness of the braid group representations in 
SU(p, q). Using just the results for 5t/(2)-theories obtained 
in Ref. 21 , infinitely many other unitary theories arise which 



have Galois conjugates satisfying Theorem IV.5 



V. COMMENTS ON NON-HERMITIAN HAMILTONIANS 

While the bulk of our paper is devoted to showing that cer- 
tain wave functions cannot be the ground states of gapped 
Hermitian Lieb-Robinson Hamiltonians, it is worth briefly 
discussing what is known about non-Hermitian Hamiltonians. 
For non-Hermitian Hamiltonians, many of the technical tools 
involving Lieb-Robinson bounds are unavailable, and so many 
results that we know in the Hermitian case are not known here. 

The first major difficulty in the non-Hermitian case is that 
even if the Hamiltonian is a sum of terms Hz which obey 
Eq. ( fT~5] >, if the Hamiltonian is not Hermitian, then we do not 
know if the Lieb-Robinson bound holds. Similarly, for a Her- 
mitian Hamiltonian, the Lieb-Robinson bound might not hold 
for evolution in imaginary time. Given that the Lieb-Robinson 
bound fails, we are also unable to prove locality of correlation 
functions in a non-Hermitian Hamiltonian even if there is a 
gap in the spectrum. 



14 



The fact that we cannot prove locality of correlation func- 
tions is relevant to the following application of the disk axiom. 
Suppose we have a Hermitian Hamiltonian which obeys the 
disk axiom with P being the projector into the ground state 
subspace. Suppose operators Ox and Oy are supported on 
small disks X and Y such that the disk axiom implies that 
POxP and POyP are both close to scalar multiples of P. 
Now, let us ask whether the operator POxOyP is also close 
to a scalar multiple of P. Consider the case in which the disks 
X and Y are far separated such that the smallest disk con- 
taining both X and Y is too large to directly apply the disk 
axiom to OxOy. Thus, the disk axiom alone does not tell 
us that OxOy is close to a scalar when projected into the 
ground state subspace. However, if we have a gapped, Her- 
mitian, Lieb-Robinson Hamiltonian then correlations decay 
exponentially in any ground state, so that POxOyP is close 
to POxPOyP, and then applying the disk axiom to POxP 
and POyP implies that POxOyP is close to a scalar mul- 
tiple of P. Unfortunately, though, in the non-Hermitian case 
we do not know that there is exponential decay of correlation 
functions, and so even if we assume the disk axiom for small 
disks, we do not see how to prove that OxOy is also close 
to a scalar when projected into the ground state subspace of a 
non-Hermitian Hamiltonian obeying the disk axiom. In fact, 
suppose we consider two states ipi , ipz such that any operator 
supported on a small disk is equal to a scalar when projected 
into the space spanned by ipi,ip2- Let us relax any require- 
ment that ipi , ip2 be ground states of a Hamiltonian, whether 
Hermitian or not, and simply take them to be arbitrary states. 
Then, we can give an example in which product operators of 
the form OxOy given above are not close to a scalar when 
projected into the ground state subspace, even though both 
Ox-,Oy are close to scalars when projected into this sub- 
space, by considering a quantum error-correcting code on a 
small number of qubits, and defining the two states ipi , ip2 on 
the large lattice by placing the qubits defining the code on far 
separated sites of the lattice, and placing all other qubits in a 
product state. 



VI. CONCLUSIONS 

To summarize, we have shown that a large class of non- 
unitary topological quantum field theories cannot be realized 
as ground states of Hermitian (quantum mechanical) Hamil- 
tonians that satisfy a Lieb-Robinson bound. This includes but 
is not limited to local and quasi-local (exponentially decay- 
ing) Hamiltonians. While our proof has been formulated for 
quantum doubles of TQFTs, it also rules out the realization of 
the constituent non-doubled TQFT in a Hermitian system: if 



the latter were to exist, it could be used to trivially construct a 
Hermitian model for the corresponding quantum double. 

The TQFTs covered by our proof include, in particular, the 
Galois conjugates of Fib and su(2)fe TQFTs for k = 3 and 
all k > 5. Among these, one case of special recent interest 
is the Yang-Lee TQFT (the Galois conjugate of Fib) underly- 
ing the proposed Gaffnian quantum Hall wave function. Our 
argument implies that this Gaffnian wave function cannot oc- 
cur as ground state of a gapped fractional quantum Hall state 
(described by a Hermitian Hamiltonian), if one considers that 
the screened Coulomb interaction satisfies a Lieb-Robinson 
bound. 
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Appendix A: Proof of Lemma[lV.3| 



We define the "quasi-adiabatic continuation operator" V s 

by 

iV s = J dtF(t)exp(iH s t)(d s H^exp(-iH s t), (Al) 

where the function F(t) is defined by (we f ollow ^, while a 
more complicated choice of F(t) was used ir l 33 l 34 j); 

— = / dMexp(-u 2 /2a 2 ), (A2) 

for t > and F(t) f or t < is defined by F(t) = -F(-t). 
The quantity a in Eq. ( |A2| > is some constant chosen below. 

We now show that d s P(s) is close to i[D s ,P(s)], and 
bound the difference between the two expressions in opera- 
tor norm. Define F(ui) to be the Fourier transform of F(t). 
One may show that 

M > AE^ \F(u))-l/w)\ < const. x(l/AE) cxp(-a 2 uj 2 /2). 

(A3) 

Let i/ji(s) denote eigenvectors of H s with eigenvalues Ei(s), 
so P(s) = E?=o M*))(M*)\- ^n, 
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n-l 



i[D s ,P(s)] = J2 E /dtF(<)|* J (s))((^(s)|exp(zi? s <)a s i/ s cxpHi/ s i)|* I ( S )))(^( S )| + / l .c. (A4) 

i=0 j>n 
n-l 

= E E F(E t -E 3 )\^( S ))((^( S )\d s H s \^(s)))(^( S )\ + h.c. 



i—0 j>n 

By linear perturbation theory, 



71—1 1 

d s p( S ) = Y / E F ^ 7 i^(s)>(<^(s)i5 5 F s |* i ( S )>)<* < ( S )i+/i.c. 



Thus, by Eqs. ( |A3|A4|A5| > 



\d s P(s) - [W s ,P(s)]\\ < const, x (\\d s H s \\/AE)exp(-a 2 AE 2 /2). 



(A5) 



(A6) 



We also have a bound on the time decay of F(t): 
\F(t)\ < const, x exp(-i 2 /2a 2 ). 



(A7) 



We now define iT>\ to be an approximation to iD s supported 
on the set of sites within distance I of X. To construct this 



approximation, we use the Lieb-Robinson bound and set 

iD\ = [ dtF(t)(d s H s ) (t), (A8) 

J\t\<l/v LR V ' 



where yd s H s j (t) is supported on the set of sites 

within distance I of X and is the approximation to 

exp(iH s t)(^d s H s ^j exp(—iH s t) given by the Lieb-Robinson 

bound. We now bound the difference ||X> S — T> l s \\. This differ- 
ence is bounded by 



dt\F{t)\\\(d s H^ (t) -exp(iff s t)(d s # s ) exp(-iH s t)\\ + j dt\F{t)\\\d s H 



'\t\<l/v LR J\t\>l/v LR 

Using the Lieb-Robinson bound and Eq. \A1\ we arrive at 

||23.-2?i|| < const.\\d s H s \\(v LR a 2 /l)(eM-l 2 /2a 2 v 2 LR ) + g(l)\X\ 



Combining Eqs. ( A6 AlOl, and choosing a = ^1/vlrAE we find that 

\\d s P{s) - \iV l s ,P(s)]\\ < const, x (\\d s H s \\/AE) [ cxp(-lAE/2v LR ) + g{l) 



(A9) 



(A10) 



(All) 



Finally, we define U(s) by U(0) = I, the identity operator, 



and d s ll(s) — iD s U(s). Since iD s is Hermitian by construc- 
tion, U(s) is unitary. Eq. ( 17 1 follows from Eq. ( Al 1 1. □ 
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